基于直接模拟蒙特卡洛方法的真空羽流不确定量化研究

doi:10.13675/j.cnki. tjjs. 190092

推进技术 ›› 2020, Vol. 41 ›› Issue (1): 73-84.DOI: 10.13675/j.cnki. tjjs. 190092

基于直接模拟蒙特卡洛方法的真空羽流不确定量化研究

陈浩¹,林震²,刘成诚¹,张斌¹,刘洪¹

1.上海交通大学航空航天学院，上海 200240;2.北京控制工程研究所，北京;100190

出版日期:2020-01-20 发布日期:2020-01-20
作者简介:陈浩，博士生，研究领域为稀薄气体动力学及高超声速气动热。E-mail：kongtianhaowen@sjtu.edu.cn
基金资助:
国家自然科学基金（NSFC-91741113；SFC-91841303）。

Uncertainty Quantification of Vacuum Plume Simulations Using Direct Simulation Monte Carlo Method

1.College of Aeronautics and Astronautics,Shanghai Jiaotong University,Shanghai 200240,China;2.Beijing Institute of Control Engineering,Beijing 100190,China

Online:2020-01-20 Published:2020-01-20

摘要/Abstract

摘要： 真实流动环境下的真空羽流必然存在着各种不确定性，那么确定性输入条件的数值模拟必然会存在偏差，因此需进一步研究不确定性对羽流流动特征的影响规律。本文采用直接模拟蒙特卡洛（DSMC）方法，对不确定性输入的羽流流场进行模拟；采用稀疏的概率配置点方法对来流、壁面及模型参数等输入不确定量进行描述，对不确定性的传播和输出目标的平均值、方差及不确定度进行计算。研究表明，流场不确定性沿流线传播至流速最大处之后迅速增强，并在声速线前出现骤减的现象；传播至声速线之后，挡板壁面输入不确定性的影响凸显。其中最为显著的是，压力不确定度在挡板驻点位置达到全场最大值，约为输入不确定量（3.54%）的2.1倍。此外，温度跳跃不确定度受到壁面温度不确定性输入的限制而近似保持为一个恒定值，约为输入不确定度的0.8倍。进而，壁面热流不确定度（5.54%）比壁面正应力不确定度（6.25%）略小，切应力不确定度最小（5.07%）。Sobol’全局敏度分析表明，喉道速度和喉道压力的输入不确定性对气动力/热不确定度的贡献是最大的，且远远超过了壁面温度和模拟分子直径不确定性输入的贡献。

关键词: 不确定量化;直接模拟蒙特卡洛法;真空羽流;配置点方法;气动特性

Abstract: Vacuum plume in real environment is uncertain whereas the numerical simulation is always deterministic, thus the influence of uncertainty on the vacuum plume is required to be elaborately studied. In this paper, the direct simulation Monte Carlo (DSMC) method was used to simulate the plume flowfield with various input uncertainties. Meanwhile, the probabilistic collocation method with sparse grid technique was used to describe the input uncertainties stemming from the incoming flow, wall condition and model parameters. The mean, variance and uncertainty of the output parameters, and the propagation of uncertainty were also calculated by the probabilistic collocation method. The corresponding results indicate that the uncertainties in flowfield sharply increase along the streamline when the flow speed reaches its maximum, and then decrease when close to the sound speed line. The influence of input uncertainty in wall temperature dominates in the downstream area of the sound speed line. Note that, the increment of uncertainty in the pressure is the most significant, which approximately equals to 2.1 times the input uncertainty (3.54%). In addition, the temperature jump uncertainty is almost constant, which is about 0.8 times the input uncertainty. This may be caused by the input uncertainty in the wall temperature that bounds the temperature jump uncertainty. As a consequence, the uncertainty in the heat flux (5.54%) is slightly smaller than that in the normal stress (6.25%), meanwhile the uncertainty in the shear stress (5.07%) is the smallest one. Finally, the Sobol’ global sensitivity analysis shows that the input velocity uncertainty and pressure uncertainty in the throat contribute the most to the uncertainty in aerodynamic characteristics, and considerably exceed the influence of the input uncertainties in the wall temperature and in the reference diameter of simulated particles.

Key words: Uncertainty quantification;DSMC;Vacuum plume;Probabilistic collocation method;Aerodynamic properties

陈浩,林震,刘成诚,张斌,刘洪. 基于直接模拟蒙特卡洛方法的真空羽流不确定量化研究 [J]. 推进技术, 2020, 41(1): 73-84.

CHEN Hao¹,LIN Zhen²,LIU Cheng-cheng¹,ZHANG Bin¹,LIU Hong¹. Uncertainty Quantification of Vacuum Plume Simulations Using Direct Simulation Monte Carlo Method[J]. Journal of Propulsion Technology, 2020, 41(1): 73-84.

参考文献

[1] 陈兵, 蔡国飙. 模拟真空羽流场的特征线法[J]. 推进技术， 2002， 23(6)： 500-504. (CHEN Bing， CAI Guo-biao. Method of Characteristics for Vacuum Plume Simulation[J]. Journal of Propulsion Technology， 2002， 23(6): 500-504.)
[2] 唐振宇，贺碧蛟，蔡国飙. 解耦N-S/DSMC方法计算推力器真空羽流的边界条件研究[J]. 推进技术， 2014， 35(7)： 897-904. (TANG Zhen-yu， HE Bi-jiao， CAI Guo-biao. Investigation on Boundary Conditions in Decoupled N-S/DSMC Method for Vacuum Plume Simulation of Thrusters[J]. Journal of Propulsion Technology， 2014， 35(7): 897-904.)
[3] 王黎珍，史纪鑫，郑世贵. 推力器真空羽流热效应计算模型修正及误差分析[J]. 航天器环境工程， 2014， 31(5)： 483-488.
[4] 贺卫东，党海燕，许明艳，等. 地外天体起飞过程真空羽流导流力热效应研究[J]. 推进技术， 2015， 36(8)： 1151-1156.
[5] 严立，王平阳，欧阳华. 月面环境发动机羽流冲击力效应模拟计算[J]. 上海交通大学学报， 2012， 46(8): 1310-1314.
[6] 蔡国飙，王慧玉，庄逢甘. 真空羽流的数值模拟及真空羽流污染的数值模拟分析研究[J]. 宇航学报， 1998， 19(3)： 1-9.
[7] 肖泽娟. 空间发动机羽流及其污染研究[D]. 上海：上海交通大学， 2007.
[8] 贺卫东. 发动机真空羽流导流力热效应的CFD/DSMC耦合仿真及试验研究[D]. 北京：北京理工大学， 2015.
[9] 张熇，蔡国飙，许映乔，等. 嫦娥三号着陆器软着陆过程中羽流仿真分析及试验研究[J]. 中国科学：技术科学， 2014， 44(4)： 344-352.
[10] 任军学，王艳，仇钎，等. 离子发动机羽流二维轴对称数值模型与验证[J]. 北京航空航天大学学报， 2011， 37(12)： 1498-1503.
[11] Crespo L G， Kenny S P， Giesy D P. The NASA Langley Multidisciplinary Uncertainty Quantification Challenge[C]. Reston， Virginia：16th AIAA Non-Deterministic Approaches Conference， 2014.
[12] Wang H， Sheen D A. Combustion Kinetic Model Uncertainty Quantification Propagation and Minimization[J]. Progress in Energy and Combustion Science， 2015， 47： 1-31.
[13] Yousefian S， Bourque G， Monaghan R D. Review of Hybrid Emissions Prediction Tools and Uncertainty Quantification Methods for Gas Turbine Combustion Systems[C]. Charlotte， North Carolina：ASME Turbo Expo 2017：Turbo machinery Technical Conference and Exposition， 2017.
[14] Beran P， Stanford B， Schrock C. Uncertainty Quantification in Aeroelasticity[J]. Annual Review of Fluid Mechanics， 2017， 49： 361-386.
[15] Bird GA. Molecular Gas Dynamics and the Direct Simulation of Gas Flows(Oxford Engineering Science Series)[M]. UK: Oxford University Press Inc.， 1994.
[16] Kulakhmetov M， Alexeenko A A. Modeling Uncertainties in Direct Simulation Monte Carlo Calculations of Hypersonic Leading-Edge Flow[J]. Journal of Spacecraft and Rockets， 2012， 49(3)： 461-472.
[17] Burt JM， Josyula E. Global Sensitivity Analysis and Uncertainty Quantification for a Hypersonic Shock Interaction Flow[J]. Journal of Thermophysics and Heat Transfer， 2015， 29(3)： 439-449.
[18] Strand J S， Goldstein D B. Global Sensitivity Analysis for DSMC Simulations of Hypersonic Shocks[J]. Journal of Computational Physics， 2013， 246： 184-206.
[19] Xiu D， Hesthaven J S. High-Order Collocation Methods for Differential Equations with Random Inputs[J]. SIAM Journal on Scientific Company， 2015， 27： 1118-1139.
[20] Bungartz H， Griebel M. Sparse Grids[J]. Acta Numerica， 2004， 13： 147-269.
[21] Lin L Y， Ren W， Zhang B. Improved Vertex Weight for Parallelization of DSMC in the Near-Continuum Regimes[J]. Journal of Aeronautics， Astronautics and Aviation， 2014， 46(2)： 88-95.
[22] Ren W， Liu H， Jin S. An Asymptotic-Preserving Monte Carlo Method for the Boltzmann Equation[J]. Journal of Computational Physics， 2014， 276： 380-404.
[23] Zhang B， Chen H， Li L Y， et al. Asymptotic-Preserving Monte Carlo Method for Large Cylinder Hypersonic Flows[J]. Journal of Thermophysics and Heat Transfer， 2018， 32(1)： 205-215.
[24] Chen H， Zhang B， Liu H. Non-Rankine-Hugoniot Shock Zone of Mach Reflection in Hypersonic Rarefied Flows[J]. Journal of Spacecraft and Rockets， 2016， 53(4)： 619-628.
[25] Zhang B， Li L Y， Mi Y X. Rarefaction Effect on the Aerodynamic Performance of Scramjet Inlet in Hypersonic Flow[J]. Journal of Aeronautics， Astronautics and Aviation， 2015， 47(4)： 357-362.
[26] 陈浩，李林颖，张斌，等. Q-K模型在氮氧离解复合反应中的评估[J]. 空气动力学学报， 2018， 36(1)： 17-21.
[27] Chen H， Li L Y， Liu C C， et al. Uncertainty Quantification of the Surface Properties in Hypersonic Rarefied Cylinder Flows[C]. Glasgow， UK：31st International Symposium on Rarefied Gas Dynamics， 2018.
[28] 刘成诚，李林颖，余彬，等. 来流速度不确定度对超声速混合层流场特性影响研究[C]. 北京：第十八届全国激波与激波管学术会议， 2018.
[29] Judd K L， Maliar L， Maliar S， et al. Smolyak Method for Solving Dynamic Economic Models: Lagrange Interpolation， Anisotropic Grid and Adaptive Domain[J]. Journal of Economic Dynamics and Control， 2014， 44： 92-123.
[30] Sudret B. Global Sensitivity Analysis Using Polynomial Chaos Expansions[J]. Reliability Engineering & System Safety， 2008， 93： 946-979.

基于直接模拟蒙特卡洛方法的真空羽流不确定量化研究

Uncertainty Quantification of Vacuum Plume Simulations Using Direct Simulation Monte Carlo Method

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 0

编辑推荐

Metrics

本文评价